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We have performed a search for muon neutrinos from dark matter annihilation in the center of 
the Sun with the 79-string configuration of the IceCube neutrino telescope. For the first time, the 
DeepCore sub-array is included in the analysis, lowering the energy threshold and extending the 
search to the austral summer. The 317 days of data collected between June 2010 and May 2011 are 
consistent with the expected background from atmospheric muons and neutrinos. Upper limits are 
therefore set on the dark matter annihilation rate, with conversions to limits on spin-dependent and 
spin-independent WIMP-proton cross-sections for WIMP masses in the range 20 - 5000 GeV. These 
are the most stringent spin-dependent WIMP-proton cross-sections limits to date above 35 GeV. 

PACS numbers: 95.35.-l-d, 14.80.Nb, 14.80.Rt, 96.50.S-, 98.70.Sa 



While the presence of dark matter (DM) in the uni- 
verse has been inferred through its gravitational inter- 
actions, its nature remains a mystery. One of the most 
promising and experimentally accessible candidates for 
DM are so-called Weakly Interacting Massive Particles 
(WIMPs) [ll, predicted in extensions of the Standard 
Model of particle physics (SM). DM may be captured 
in large celestial bodies like the Sun [2| where self- 
annihilation to SM particles can result in a flux of high- 
energy neutrinos. These neutrinos can be searched for 
as a point-hke source by IceCube 1,11. These indirect 
searches for DM are sensitive to the WIMP-proton scat- 
tering cross section, which initiates the capture process 
in the Sun. They complement direct DM searches on 



Earth as they scale with the averaged DM density along 
the solar circle and are more sensitive to low WIMP ve- 
locities 0. Indirect searches depend onlyweakly on the 
underlying WIMP velocity distribution [6| and we have 
chosen parameters to be conservative in our analysis. 

In this work, wc present new IceCube limits on dark 
matter captured by the Sun, with data taken in the 79- 
string configuration of the detector. This analysis incor- 
porates two significant additions compared to previous 
work. Firstly, we extend the search to the southern hemi- 
sphere. This doubles the livetime of the analysis, but 
imposes new challenges to reduce the downgoing atmo- 
spheric muon background. Secondly, we search for neu- 
trinos from WIMPs with masses (m^) as low as 20 GeV 
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whereas past IceCube searches have only been sensitive 
from 50 GeV. 

The IceCube detector 0] is situated at the South Pole. 
Digital Optical Modules (DOMs) arranged on vertical 
strings deep in the ice sheet record the Cherenkov 
light emitted by relativistic charged particles created 
in neutrino interactions in the ice. The detection of 
photon yields and arrival times in DOMs allows the 
reconstruction of direction and energy of the secondaries. 
This analysis used 317 live-days of data taken between 
June 2010 and May 2011. During this period, the 
detector was operating in its 79-string configuration, 
which includes six more densely instrumented strings in 
the center of the array, optimized for low energies. These 
strings feature reduced vertical spacing between DOMs 
and higher quantum efficiency photomultiplier tubes. 
Along with the seven surrounding regular strings, they 
form the DeepCore subarray Both the improvement 
in livetime and in energy threshold, which this analysis 
has achieved over previous IceCube analyses, can be 
attributed to the use of the DeepCore array. 

All signal simulations are made with DarkSUSY Q and 
WimpSim [lo| . which describe the capture and annihila- 
tion of WIMPs inside the Sun and the consequent pro- 
duction, interaction, and propagation of neutrinos from 
the core of the Sun to the detector, including three-flavor 
oscillations and matter effects. The primary WIMP an- 
nihilation spectrum is very model dependent, owing to 
different branching ratios into SM particles. We approx- 
imately bracket the range of possible models by assum- 
ing 100% branching into two channels with very differ- 
ent characteristics. The 'hard' channel, XX ~^ W~^W~ 
(r+r" below m^ — 80.4 GeV) and the 'soft' channel 
XX ~^ For this work, WIMP masses ranging from 
20 GeV to 5 TeV arc simulated. 

The background in this search consists of muons 
and neutrinos created in cosmic ray interactions in the 
Earth's atmosphere. The dominant down-going muon 
component is simulated with CORSIKA [llj, including sim- 
ulations of single and coincident air showers. The and 
i/g components of the atmospheric spectrum are gener- 
ated following the Honda flux model 12] . For verification 
and cross-checks, a dedicated simulation of atmospheric 
us below 200 GeV is performed with GENIE [l^ . The total 
background at final analysis level from solar atmospheric 
neutrinos, originating from cosmic ray interactions in the 
Sun's atmosphere, has been calculated to be of order 1 
event, depending on the flux model 14-1^. To reduce 
the dependence on simulation and associated systematic 
errors, we use off-source data to estimate the background 
at all analysis levels. 

Propagation of muons through the ice is simulated [l^l , 
and transport of light from these particles to the DOMs 
is performed using direct photon tracking [lij . taking 
into account measured ice properties (20j . Particle and 



photon propagation simulations at the lowest targeted 
energies below 50 GeV have been independently verified 
using GEANT4 

In this work the full dataset is split into three indepen- 
dent non-overlapping event selections; first into 'summer' 
and 'winter' seasons, when the Sun is above and below 
the horizon, respectively. The 'winter' dataset is further 
split into a low energy sample (WL), with focus on start- 
ing neutrino-induced muon tracks inside DeepCore, and 
a higher energy sample (WH), aiming to select track- 
like events with no particular containment requirement. 
The 'summer' selection is a dedicated low energy event 
sample (SL) for which the surrounding IceCube strings 
are used as an active muon veto in order to select start- 
ing neutrino-induced events within DeepCore. Separa- 
tion into these samples is necessary owing to the differ- 
ent characteristics of the overwhelming down-going muon 
background within each dataset. The event selection is 
carried out separately for each independent sample and 
the final search is conducted using a combined likelihood 
function. In order to avoid potential bias, a strict blind- 
ness criterion is imposed by scrambling the azimuthal 
position of the Sun in data. 

In IceCube, filters pre-select data to enhance the con- 
tent of signal-like muon events above the dominant at- 
mospheric muon background. To increase the signal-to- 
background ratio, we only select events that pass any 
of three filters: the dedicated DeepCore low energy fil- 
ter and two filters selecting muon-like events with an 
upwards pointing track reconstruction. At this point, 
the dataset is split into the two seasonal streams, where 
September 22nd 2010 and March 22nd 2011 mark the 
beginning and ending dates of the SL selection. We first 
discuss the additional 'winter' cut selections: Cuts are 
applied on the zenith angle and quality of the likelihood- 
based track reconstruction, on hit and string multiplic- 
ity, and on timing and topological variables. For Deep- 
Core contained events, the zenith acceptance region is 
extended to reflect the broadened signal point spread 
function at low energies. 

The first data reduction is followed by additional pro- 
cessing, including an estimate of the angular uncertainty 
of the muon track fit. Some signal neutrinos will ar- 
rive in coincidence with atmospheric backgrounds. In 
order to retain these signal events, a set of topological 
criteria are applied to 'split' these combined hit patterns 
into distinct sub-events. These sub-events are then pro- 
cessed as above, and undergo all subsequent event selec- 
tion in their own right. Following the addition of events 
from splitting, the dataset is divided into independent 
low and high energy event samples. For events to be in- 
cluded in the WL sample, we demand that the number 
of hit DOMs within DeepCore must be larger than out- 
side. Additionally, the number of outside hits must be 
less than seven. This ensures that events with a long lever 
arm and therefore good angular resolution are assigned 



4 




0.99 



0.992 



0.994 



0.996 



0.998 1 



FIG. 1. Cosine of the angle between the reconstructed track 
and the direction of the Sun, vp, for observed events (squares) 
with one standard deviation error bars, and the atmospheric 
background expectation from atmospheric muons and neutri- 
nos (dashed line). Also shown is a simulated signal (1 TeV 
hard for dataset WH, 50 GeV hard for datasets WL and SL) 
scaled to nf^ (details in Table 



to the complement sample. Events that fail the WL eri- 
teria are classified as WH events, and undergo a series of 
additional, stricter cuts on the same variables as in the 
initial event selection. WL events, conversely, undergo a 
veto cut, removing events with hits in the 10 upper most 
layers of DOMs on the regular (non-DeepCore) strings. 
The final background reduction utilizes one Boosted De- 
cision Tree (BDT) [39] for each dataset to discriminate 
true up-going muon-like events from mis-reconstructed 
atmospheric muons. For training and testing an inde- 
pendent, high statistics, set of signal simulations is used 
and discarded afterwards. For background training this 
multivariate analysis uses one month of data recorded 
when the Sun was outside the analysis region. Through 
an iterative process, individual variables were removed 
and added and the performance of the BDT evaluated, 
until we arrived at a final set of 14 variables in the WH 
stream and 10 in the WL stream. All input distributions 
for simulated backgrounds and data are in good agree- 
ment. The selected variables describe both the quality 
of the track reconstruction and the time evolution of the 
pattern of hit DOMs and spatial positions within the de- 
tector. 

The SL event sample uses a different set of cuts, be- 
cause the dominant background is comprised of well- 
reconstructed down-going muons penetrating the detec- 
tor. To reduce these backgrounds, we focus on low 
signals with a reconstructed neutrino interaction vertex 



inside the DeepCore fiducial volume. Selecting only these 
events, cuts are placed on the zenith angle of the track 
reconstruction, hit multiplicity, and vertical extension of 
the event. A 14 DOM layer top veto is imposed to re- 
ject down-going events. Additionally, events are required 
to be DeepCore dominated (defined in the same way 
as for the 'winter' analysis) and fulfill a tight hit-time 
containment criterion. The final step in background re- 
jection again consists of one BDT with 10 input vari- 
ables. These are selected using the same iterative selec- 
tion process. Track quality parameters yield less separa- 
tion power within this down-going sample. As a result, 
the final BDT input observables mainly describe the de- 
gree of containment and the vertical and lateral extension 
of the event within the detector. 

The cut on the BDT score is optimized for each 
event selection to minimize the model rejection factor, 
MRF [2^, in the full likelihood analysis. The final step 
of the analysis is a likelihood ratio hypothesis test based 
on the values of the reconstructed angle to the Sun 
using the Feldman-Cousins unified approach This 
results in confidence intervals for the mean number of 
signal events, fis- The required probability densities for 
signal are computed from simulations, while for back- 
ground they are based on real data events at the final 
selection level, with scrambled azimuth direction. A sin- 
gle result is calculated from all three data samples with 
a combined likelihood, constructed from the set of three 
independent probability distributions of signal and back- 
ground, weighting each by the respective livetime and 
effective volume (see Ref. [J| for details). 

After unblinding the direction of the events in the final 
data samples, the observed distributions are compared to 
the expected background distributions from atmospheric 
muons and neutrinos, shown in Fig. [TJ The observed 
number of events from the direction of the Sun are con- 
sistent with the background-only hypothesis. Upper 90% 
CL limits on /Xs are calculated and listed for each signal 
hypothesis in Table HI 

The upper limit on Hs can be translated into a limit on 
the signal fiux and annihilation rate in the Sun. The ef- 
fect of different sources of systematic uncertainties on sig- 
nal flux expectations is calculated for three signal energy 
regions, defined in Table HIl bv corresponding benchmark 
WIMP masses. Sources of uncertainties are divided into 
two classes; measurement and parameterization errors on 
cross sections and neutrino properties on the one hand 
and limitations in the detector simulation and uncertain- 
ties in detector calibrations on the other hand. The first 
class, Class-I, affects signal normalizations only, whereas 
the latter (Class-II) alters signal acceptance and intro- 
duces changes in the point spread function that is the ba- 
sis for the likelihood analysis. Class-II uncertainties are 
evaluated using alternative signal simulations with varied 
calibration parameters, processed through the same anal- 
ysis chain, and evaluated with the full multi-dataset com- 
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TABLE I. Results from the combination of the three independent datasets. Upper 90% hmits on the number of signal events 
/Xg", the WIMP annihilation rate in the Sun Fa, the muon flux and neutrino flux and the WIMP-proton scattering 
cross-sections (spin- independent, (Jsi,p, and spin-dependent, asD.p), at the 90% confidence level including systematic errors. 
The sensitivity "1>^ (see text) is shown for comparison. 
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bined likelihood. This procedure explicitly determines 
the systematic effect on /i^. 

Uncertainties in neutrino-nucleon cross-sections arise 
in the parameterization of the CTEQ6-DIS parton dis- 
tribution functions as used in nusigma [lH. In addi- 
tion to this theoretical uncertainty on a^^ the energy de- 
pendent error on the experimental (Ti^-measurement (23j 
is included. The uncertainty in neutrino oscillation 
parameters used in signal flux calculations is investi- 
gated through variations of mixing parameters within the 
quoted la regions Here, the dominant effect results 
from the least constrained mixing angle, 623, maximizing 
tau (dis)appearance within the expected flux expecta- 
tion. 

The second class of uncertainties includes absolute 
calibration and DOM to DOM variation of sensitivity, 
optical properties of the glacial ice, and photon prop- 
agation to the detector. The systematic uncertainties 
on absolute DOM sensitivity are evaluated with sets of 
signal simulations with an overall shift of 10% in DOM 
efficiency. As baseline simulations do not account for 
varying relative DOM efficiency, dedicated signal simu- 
lations are performed with individual DOM efficiencies 
from a Gaussian fitted to the in-situ measured spread 
(ct = 0.087) and centered around the nominal value. 
Optical properties of the glacial ice are measured ^ 
and characterized in models that are parameterizations 
of the absorption and scattering coefficients 0jS 0j func- 
tion of depth and position in the detector. Two such 
models 20|, |21|, differing in parameterization techniques, 
are considered to bracket the uncertainty in light yield 



resulting from the ice description. Individual uncertain- 
ties, listed in Table HIl are added in quadrature to obtain 
the total systematic uncertainty for each benchmark 
mass region. 



The upper limits on Hs for each signal hypothesis are 
then converted to limits on the neutrino to muon con- 
version rate and, through DarkSUSY to limits on 
the WIMP annihilation rate in the Sun, Pa- For bet- 
ter comparison to other experiments limits on the neu- 
trino flux (^i/) from the Sun, and the corresponding in- 
duced muon flux in the ice ($/i), both integrated above 
1 GeV, are computed at the 90% confldence level. These 
limits are listed in Table HI Also specified is the me- 
dian sensitivity, 4>^, derived from simulations without 
signal. Under the assumption of equilibrium between 
WIMP capture and annihilation in the Sun, limits on Pa 
are converted into limits on the spin-dependent, crgu p, 
and spin-independent, asi^p, WIMP-proton scattering 
cross-sections, using the method from Ref 0. Results 
are listed in Table |T] and shown in Fig. [2] together with 
other experimental limits 28-31]. We assume a stan- 
dard DM halo with a local density of O.SGeV/cm^ ^ 
and a Maxwellian WIMP velocity distribution with an 
RMS velocity of 270km/s. We do not include the de- 
tailed effects of diffusion and planets upon the capture 
rate, as the simple free-space approximation Q included 



in DarkSUSY is found to be accurate l41|. Limits on the 
WIMP-nucleon scattering cross section can also be de- 
duced from limits on mono-jet and mono-photon signals 
at hadron colliders, but these depend strongly on the 
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I 1 MSSM incl. XENON (2012) ATLAS + CMS (2012) 

' : DAMA no channeling (2008) 

CDMS(2010) - 

CDMS2keV reanalyzed (2011) 

CoGENT(2010) 

XENON100(2012) _ 




lceCube2012(bB) 
lceCube2012(W'W")* 
* (tV tor < m^^ = 80.4GeV) 

^ \ I I 



TABLE 11. Systematic errors on signal flux expectations in 
percent. Class-II uncertainties marked * 



logic ( m /GeV) 



] MSSM incl. XENON (2012) ATLAS + CMS (2012) 
DAMA no clianneling (2008) 

■ COUPP(2011) 

■ KIMS (2012) 

- Simple (2011) 

■ P[CASSO(2012) 

■ SUPER-K(2011)(bB) 

- SUPER-K(2011)(WV) 



-V. 



lceCube2012(bB) 
— lceCube2012(W^'\/y 
» for m < = 80.4GeV) 

, f , , L 




Iog10 ( m / GeV) 



FIG. 2. 90% CL upper limits on crsi,p (top figure) and 
o'SD.p (bottom figure) for iiard and soft anniliilation ciian- 
nels over a range of WIMP masses. Systematic uncertain- 
ties are included. The shaded region represents an allowed 
IVISSIVI parameter space (MSSM-25 [H]) taking into account 
recent accelerator [231, cosmological and direct DM search 
constraints. Results from Super-K [2I, COUPP PI- 
CASSO [31II KIMS [13, CDMS [32,l3lL XENON- 100 [13, 
CoGeNT [3i], Simple [si] and DAMA [si. Hi] are shown for 
comparison. 



choice of the underlying effective theory and mediator 
43-44 1, and consequently not included in Fig. [2] 



masses 



In conclusion, we have presented the most stringent 
limits to date on the spin-dependent WIMP-proton cross- 
section for WIMP masses above 35 GeV. With this 
dataset, we have demonstrated for the first time the abil- 
ity of IceCube to probe very low mass dark matter mod- 
els. This has been accomplished through effective use 
of the DeepCore sub-array. Furthermore, we have ac- 
cessed the southern sky for the first time by incorpo- 
rating strong vetos against the large atmospheric muon 
backgrounds. The added livetime has been shown to im- 
prove the presented limits. IceCube has now achieved 
limits that strongly constrain dark matter models and 
that will impact global fits of the allowed dark matter 



Source 


mass 


ranges 


(GeV) 




< 35 


35-100 


> 100 


v oscillations 


D 


D 


D 


j/-nucleon cross-section 


7 


5.5 


3.5 


/^-propagation in ice 


<1 


<1 


<1 


Time, position calibration 


5 


5 


5 


DOM sensitivity spread* 


6 


3 


10 


Photon propagation in ice* 


15 


10 


5 


Absolute DOM efficiency* 


50 


20 


15 


Total uncertainty 


54 


25 


21 



parameter space. This impact will only increase in the 
future, as analysis techniques improve and detector live- 
time increases. 
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